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ABSTRACT 

We compute the redshift space power spectrum of two X-ray cluster samples: 
the X-ray Brightest Abell Cluster Sample (XBACS) and the Brightest Cluster Sam- 
ple (BCS) using the method developed by Feldman, Kaiser & Peacock. The power 
spectrums derived for these samples are in agreement with determinations of other 
optical and X-ray cluster samples. For XBACS we find the largest power spectrum 
amplitude expected given the high richness of this sample (i? > 2). In the range 
0.05/iMpc~^ < k < 0.4/iMpc~^ the power spectrum shows a power law behavior 
P{k) oc fc" with an index n ~ —1.2. In a similar range 0.04/iMpc^^ < k < 0.3/iMpc^^ 
BCS power spectrum has a smaller amplitude with index n ~ —1.0. We do not find 
significant evidence for a peak at /c ~ 0.05/iMpc~^ suggesting that claims such of fea- 
ture detections in some cluster samples could relay on artificial inhomogeneities of the 
data. We compare our results with power spectrum predictions derived by Moscardini 
et al. within current cosmological models (LCDM and OCDM). For XBACS we find 
that both models underestimate the amplitude of the power spectrum but for BCS 
there is reasonably good agreement at A; ^ 0.03ft-Mpc~^ for both models. 

Key words: galaxies: clusters: general-large scale structure of Universe-X-rays: galax- 
ies 



1 INTRODUCTION 

The distribution of matter at very large scales can be traced 
using galaxy clusters that are the largest virialized objects 
in the Universe. This distribution is deeply connected to 
the fluctuations in the primordial density field since at very 
large scales gravitational effects are still linear. Assuming a 
Gaussian distribution of fiuctuations the two point correla- 
tion function ^(r) or the power spectrum P{k) are statistical 
tools suitable to give a complete description of the matter 
distribution. From a mathematical point of view any of these 
functions are equivalent since they form a Fourier transform 
pair. In the last years, considerably efi'ort has been carried 
out applying these statistics to difi'erent observational sam- 
ples of galaxy clusters. 

The Abell (1958) catalogue and its extension Abell- 
ACO (Abell, Corwin & Olowin, 1989), constructed by vi- 
sual inspection of Palomar photographic plates, is prob- 
ably the most widely used cluster survey. In pioneering 
works Peacock & West (1992) and Jing & Valdarnini (1993) 
have computed the power spectrum for different Abell clus- 
ter samples. At large wave-numbers k ^ 0.05^Mpc~^ they 



found a power law behavior P{k) oc k" with an index 
— 1.4 <n<. —1.7. At lower wave-numbers, k < 0.05/iMpc~^, 
there is a striking flattening of the power spectrum chang- 
ing dramatically the value of n. More recently, Einasto et 
al. (1997) and Retzlaff et al. (1998) using other AbeU-ACO 
samples found a peak at fc ^ 0.04 — 0.05/iMpc~^ near the 
wave-number where the flattening starts. Miller & Batuski 
(2000) reanalyzing an AbeU-ACO sample claim that this 
peak is only present due to the inclusion of clusters of rich- 
ness R = and clusters with estimated (not measured) red- 
shifts. They also point out that R = clusters are outside 
the statistical sample of Abell (1958) and do not constitute 
a fair sample. It should be noted the importance of such 
features due to the deep link with primordial mass power 
spectrum of current cosmological models. 

It has been claimed that projection effects can be 
strongly present in optical cluster catalogue and could in- 
crease artificially their power spectrum amplitude. One at- 
tempt to minimize subjectivites present in optically selected 
cluster catalogue is the identification of clusters in pho- 
tografic plates by means of automatic computer algorithms 
like in the APM cluster survey (Dahon et al. 1997). The 
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power spectrum obtained by Tadros, Efstathiou & Dalton 
(1998) for this cluster sample shows a smooth transition 
from positive to negative slopes near k = 0.03 — 0.04/iMpc~^ 
without the detection of any peak. Another way to avoid this 
problem is the identification of clusters in the X-ray band 
where the highly peaked emission in cluster centers mini- 
mizes the probability of projection effects. Bohringer et al. 
(1998) present the ROSAT ESO Flux Limited X-ray (RE- 
FLEX) survey that is an optically confirmed X-ray cluster 
sample in the southern hemisphere. For a relatively small 
subsample they also present one of the highest power spec- 
trum amplitude for a cluster sample without detecting the 
presence of a peak. There are another two X-ray flux lim- 
ited cluster catalogue constructed from the ROSAT all-sky 
survey (RASS): the X-ray Brightest Abell-type Cluster sam- 
ple (XBACS; Ebehng et al. 1996) and the Brightest Cluster 
Sample (BCS; Ebeling et al.l998). Recently, an extension 
of BCS to a lower flux limit has been released forming the 
extended BCS (EBCS, Ebeling et al. 2000). This sample 
is constructed from the extended X-ray emission of RASS, 
while the XBACS is based on previously optically selected 
clusters. In a recent paper, Moscardini et al. (2000) pro- 
vide predictions for the correlation function and power spec- 
trum of X-ray flux limited surveys for different cosmological 
models. They apply their models to XBACS, BCS and RE- 
FLEX following the non-linear evolution of clustering and 
using theoretical and observational relations between intrin- 
sic properties of clusters like mass, luminosity and temper- 
ature. 

The aim of this paper is to estimate the redshift space 
power spectrum of XBACS and BCS in order to compare 
with observational results for different samples and model 
predictions. The paper is organized in the following way: 
In section 2 we describe the X-ray cluster samples to be 
analyzed. In section 3 we present the method applied to 
estimate the power spectrum of a flux limited sample with 
the derivation of the corresponding errors. In section 4 the 
main results are discussed and we present the conclusions in 
section 5. 



2 XBACS AND BCS 

In this section we briefly describe the main characteristics 
of the samples relevant for the power spectrum estimation. 
The XBACS comprise 242 Abell-ACO X-ray confirmed clus- 
ters distributed in the whole celestial sphere excluding the 
strip of low galactic latitude \biim\ < 20°. The flux limit of 
the sample is Fum — 5 x 10"^^ erg s~^ cm~^ in the en- 
ergy range 0.1-2.4 keV and the redshift limit is z = 0.2 . 
Ebeling et al. (1996) estimate that the overall completeness 
of the XBACS is ^ 80% at the above Fum- This sample is 
free of volume incompleteness since is constructed from the 
brightest Abell/ACO clusters which suffer of incompleteness 
problems only for poorest clusters (Plionis & Kolokotronis 
1998). 

The original BCS consists of 201 clusters distributed in 
the northern equatorial hemisphere {5 > 0°) with the same 
cut in galactic latitude \biim\ < 20° and a redshift limit 
z < 0.3. The corresponding flux limit is similar to XBACS 
{Fiim = 4.4 X 10~^^erg cm~^) in the same band with 
a completeness of ~ 90%. The extension of BCS (Ebeling 
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Figure 1. Altoff projection using galactic coordinate of the dis- 
tribution of XBACS (crosses) and BCS (original BCS are shown 

as open circles and its extension to a lower flux limit as filled 
squares). The dashed lines show the limit |6Kml ^ 20° and the 
dotted line 5 = 0° 

et al. 2000) includes 100 additional clusters up to a lower 
limit Fiim = 2.8 x 10~^^ erg s~^ cm~^. These two samples 
form jointly the EBCS comprising 301 clusters with a nomi- 
nal completeness relatively low ^ 75%. Ebeling et al. (1996, 
1998 and 2000) give right ascensions, declinations (J2000.0) 
and redshifts for XBACS, BCS and its extension, respec- 
tively. 

In Figure 1 we show the angular distribution in galactic 
coordinates of XBACS (crosses) and BCS (open circles for 
the original sample and filled squares for its extension). The 
superposition of open circles and crosses indicate that > 50% 
of BCS are included in the XBACS. 



3 POWER SPECTRUM ESTIMATION 
3.1 Introduction 

In this section we outline the scheme applied to estimate 

the redshift space power spectrum. If a redshift survey is 
used to compute the power spectrum of a distribution, then 
the resulting power spectrum is the convolution of the real 
spectrum with the catalogue window function. This window 
function depends on the geometry of the survey and on the 
technique employed to calculate it and should be as narrow 
as possible. Several authors present different methods to es- 
timate the redshift space power spectrum of a sample. In a 
pioneer work of Feldman, Kaiser & Peacock (1994, hereafter 
FKP) provide a suitable method for a flux limited sample 
with a detailed analysis of errors. Tegniark (1995) present 
another way to calculate the power spectrum that maxi- 
mizes the spectral resolution for different survey geometries. 
Tadros & Efstathiou (1996) give a variation of FKP, but for 
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BCS 


3.13 


1.59 


16.22 


74.83 
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a volume limited sample. Analyzing different methods to 

compute the power spectrum, Tcgmark et al. (1998) show 
that FKP method is the most appropriate for wave-numbers 
greater than the point where the value of the window func- 
tion decrease a factor 2 of its value at the origin. We apply 
the method derived by FKP in the version described by 
Hoyle et al. (1999). This method is suitable for our samples 
due to its simple quasi-spherical geometry (Sutherland et al. 
1999). 



3.2 Spatial Distribution 

We have used the standard transformation from redshift z to 
comovirig distance r for a model universe with a dimension- 
less density parameter equal to unity and a null cosmological 
constant (Mattig 1958): 

r = 2-^(1 -(l + z)-i/2) (1) 

where c is the speed of light and Ho = 100 h km s~^/Mpc 
is the Hubble constant. 

In the upper panel of Figure 2 we show the histogram 
of the distribution of XBACS (solid line) and BCS (dashed 
line) as a function of comoving distance. We have taken bins 
of comovirig distance width dr ~ 35/i^^Mpc. The smooth 
curves are the fitting formulae proposed by FKP: 



N{r) 



^ max ( ~ 



[l + (- 



(2) 



with parameters quoted in Table 1 by a chi-square maximum 
like-hood method. 

An often applied method to test for systomatics gradi- 
ents in the samples is to compute the observed comoving 
number density 



n(r) 



N{r) 



(3) 



where N{r) is the observed number of clusters in a bin 
at a distance r and Q is the solid angle covered by the sam- 
ple. Actually, the sky coverage is a function of the flux F, but 
this function is not available for these catalogues. Following 
Moscardini et al. (2000) we assume that the actual sky cover- 
age is a constant, for XBACS is = 47r(l — sirz(6iim)) = 8.27 
and for BCS is ^ = 4.13. 

Since these are a flux limited samples we have normal- 
ized n{r) to the expected comoving number density no{r) 
obtained integrating the luminosity function. 



no(r) 



/: 



AL- 



(4) 



where Lmin ~ 47rr^(l -|- z)Fii^. We have adopted the 
Schechter luminosity function parameters from Plionis & 
Kolokotronis (1998) and we quote them in Table 2. 




Figure 2. Upper panel: Histogram of the spatial distribution of 
XBACS (solid line) and BCS (dashed line) . The smooth curves are 
the corresponding fitting formulae given by equation (2). Lower 
panel: Normalized comoving number density of XBACS (filled 
circles) and BCS (open circles) along the radial direction. 



Table 2. Parameters of Schechter luminosity function. A 
units are W-^h^MTpc-^{lQ'^'^h-'^ev^-'^)"-'^ . L* units are 
1044/i-2ergs-i 



sample 


A 


L* 


a 


XBACS 


1.955 


1.048 


1.21 


BCS 


1.246 


2.275 


1.85 



In the lower panel of Figure 2, we plot the normalized 

comovirig number density n(r)/no{r) of XBACS (filled cir- 
cles) and BCS (open circles). This plot shows the absence of 
spatial gradients in the samples. See Plionis & Kolokotronis 
(1998) for a detailed description of these effects. 



3.3 The method 

If the cluster sample has Nc clusters with vector positions 
Xc we construct a random catalog with the same geometry 
and selection function of the sample with Nr points with 
vector positions x,.. Then, we define a quantity with mean 
value zero as: 

<5(k) = D(k) - aWOs.) (5) 
where 

D(k) =^a;(a:c)e'''"'= (6) 

c=l 

is the Fourier transform of the cluster distribution, and 
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(7) 



is the Fourier transform of the window function of the sur- 
vey. In Equation (5) a = Sc/Sr with 



Sc = [xc] and Sr 



(8) 



Assuming Gaussian density fluctuations, FKP derive a 
weight function 

1 



uj{r) 



l + n(r)P^(fc) 



(9) 



that minimizes the power spectrum variance. To compute 
these weights the actual value of the power spectrum is 
needed. To solve this problem we propose different values 
for Pw{k) as an initial guess. Having defined 5(k) by equa- 
tion (5) the power spectrum estimator is obtained by: 



P{k) = {\S{k)f -a{l + a)Sr)/C 

where C is a normalization constant defined by: 



(10) 




. P(k) extracted from the Hubble volume 
P(k) average from 27 mocks (Hubble volume)(FKP) 



k/{Mpc-i h) 



(11) 



and V is the volume where periodicity is assumed. We 
have adopted the above definition of a = Sc/Sr in or- 
der to recover the definition of P{k) given by equation 
(2.4.5) of FKP. Finally, assuming isotropy we compute the 

power spectrum estimator averaging over spherical shells 
k < < k + dk where there are iV^ wavenumber vectors 
ki: 



P{k) 



1 



iVfc 



^P(k. 



(12) 



We compute spectral densities at the nniltiples of the 
fundamental mode in order to avoid oversampling of the 
spectra that could produce spurious features. In order to 
compute the Fourier transform of the distribution of the 
cluster sample (equation (6)) and of the random catalogue 
(equation (7)) we use a Fast Fourier Transform (FFT) al- 
gorithm (Press et al., 1986). We compute these quantities 
embedding the distributions within a periodic larger cubic 
volume V = r^^^ divided in A'^ cells per side. To use FFT 
we assign the spatial distribution of points (clusters or ran- 
dom) into the grid by means of different weight assignment 
schemes. 



3.4 Error estimations 

We have estimated errors of the power spectrum using equa- 
tion 2.4.6 of FKP: 

= il? I^WQ(ki - k,) - S(k, - k,)f (13) 

i=l 3 = 1 

where 



Q(k) = a^n(a;^)w^(a;^)e*-"7C 



(14) 



Figure 3. The solid line shows the power spectrum of haloes 
identified in the complete Hubble Volume Simulations and points 

arc the average haloes power spectrum of 27 mocks catalogues ex- 
tracted from tlic same simulation. The later was computed using 
the FKP formalism (eq.lO). The error bars show the correspond- 
ing l-<7 standard deviation. 



and 



5(k) = a(l + a)J2 w^(a;r-)e*'' "7C. 



(15) 



In equation (13) ki and kj are assumed to belong to 

the same spherical shell < |k| < fe 4- dk. Note that ^(k — 
(0,0,0)) = a(l -|- a)Sr is the second term of the definition 
of P{k), equation (11). 



3.5 N-body Simulations 

In the previous section we have used equation 2.4.6 of FKP 
in order to estimate power spectrum errors. To test the er- 
ror analysis of FKP we have used a second method based 
on the variance of P{k) obtained from mock XBACS cata- 
logue. We have constructed these mock catalogue from the 
Hubble volume simulations carried out by the Virgo Consor- 
tium (Jenkins et al. 1998). The sirimlation use lO''^ particles 
in a cubic volume of 3000 h^^ Mpc per side resulting in a 
mass per particle of 2.25 x 10^^h~^MQ . The initial condi- 
tion were generated using the Lambda Cold Dark Matter 
(LCDM) model with the following parameters: dimension- 
less matter density Sim = 0.3, dimensionless cosmological 
constant density Qa = 0.7, h = 0.7 and relative mass fluctu- 
ation in a sphere of 8 Mpc as = 0.9. This normalization 
is consistent both with the observed abundance at z = 
of rich clusters and the fluctuations detected by the satel- 
lite COBE. From this simulation, we have used the haloes 
of particles kindly provided by Caxlton Baugh identified us- 
ing a standard "friends of friends" algorithm. The selected 
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Figure 4. Dependence of XBACS (upper curves) and BCS (lower 
curves) window functions on different parameters (see labels). 
Panel A: weight assignment scheme. Panel B: number of grids A''. 
Panel C: box side ri,ox Panel D: guessed power spectrum P^{k). 
To maJje the plot clearer wc have lowered the amplitude of BCS 
curves by one order of magnitude. 



Figure 5. The estimated power spectrum of XBACS (upper 
curves) and BCS (lower curves). The different panels correspond 
to different mass assignment scheme, number of grids, weights 
and box size as indicated in Figure 3. The logarithmic amplitude 
of BCS power spectrum is shifted -0.5 to make the plot clearer. 



haloes have a mean separation of dc = 50/i~^ Mpc in order 
to reproduce the same number of cluster of the catalogue. 
This choice is also consistent with the correlation length 
of ro ~ 21/i"^ (Abadi, Lambas & Muriel, 1998) from the 
ro — dc relation obtained by Colberg et al. 2000. The volume 
of the simulation allows to extract 27 (3 per side) indepen- 
dent boxes with a side comparable to the XBACS cluster 
catalogue. In order to mimic some of the observational con- 
strains, we have applied to each one of these boxes the same 
cuts in redshift and galactic latitude of the XBACS sample. 
Also, we convert center of mass positions of haloes from real 
to redshift space using the corresponding peculiar velocity 
of the halo (the conversion from comoving distance to red- 
shift and viceversa was realized using equation 1). Finally, 
we convert total masses to X-ray luminosities using the em- 
pirical relation M — Lx proposed by Reiprich & Bohringer 
2000: 



M 



= 4.7x 10'*(: 



(16) 



We remove clusters with X-ray luminosity lesser than 
the luminosity limit which is given by (Peacock 2000) 

47r 



-dL(z)Fii. 



erg s ^cm ^ in the energy range 0.1- 



1+2 

where Fum ~ 5x10^^^ 
2.4 kcV and dh is the luminosity distance 

dL = ril + z) 



(17) 



(18) 

with r given by equation (1). 

Using the method described in section 3 we compute 
the power spectrum of each one of the 27 constructed mock 



catalogue and then we compute the mean power spectrum 
and the corresponding la standard deviations. Using this 
model and parameters the haloes obtained have both a red- 
shift space distribution and a power spectrum similar to 
XBACS. 

In order to test the ability of the method to reproduce 
the correct power spectrum of clusters we compare the power 
spectrum of the haloes identified in the Hubble Volume Sim- 
ulation with the average power spectrum of the 27 mocks 
catalogues using equation 10. We have computed the power 
spectrum of the Hubble Volume haloes using 



P{k) 



IE' 



1 



(19) 



where Vh is the cubic volume of the Hubble simulation and 
Nfi is the number of haloes identified in this volume. 

In figure 3 the solid line is the power spectrum of Hub- 
ble Volume haloes and the points is the average haloes power 
spectrum of the 27 mocks catalogues. The error bars show 
the corresponding l-cr standard deviations. For this compar- 
ison we have chosen haloes identified using dc = 30Mpch~^ 
in order to obtain a large number of haloes in each mock 
(~ 1000). This figure shows that the agreement between 
both estimates is very good for k > 0.05hMpc~^ but for 
larger scales the power spectrum amplitude of the mocks 
catalogues is slightly lower. The error bars in the all range 
make the determinations undistinguishables showing that 
the method derived from FKP is a good estimator of power 
spectrum of a spatial distribution. 

We also note that the method can take correctly into 
account the galactic extinction which can generate an over- 
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estimate of the amplitude in the power spectrum mainly on 
large scales (Vogeley 1998). 



4 DISCUSSION 

First, we compute the window function of XBACS and BCS 
using equation (7). In Figure 4 wc plot the window func- 
tion module square {\W{k)\/Nr)'^ normalized to the number 
of points used to generate the random catalogue. For both 
samples, these random catalogue have A^^ = 10® points dis- 
tributed within the same geometrical limits of each sample 
and with the radial number density according to equation 
(3). In order to avoid superposition of XBACS and BCS 
curves in this Figure the amplitude of {W{\k\)/Nr)'^ for BCS 
is plotted with a factor 0.1. The window function can be very 
well approximated by a power law with index n ~ —4 over 
the range O.OlftMpc^^ < k < O.l^Mpc^^ for both samples. 
This steep negative slope guarantees the validity of equa- 
tion (10) (see Tadros & Esfthatiou 1996). In the first three 
panels we test different parameters involved in the appli- 
cation of the FFT technique: panel A, weight assignment 
scheme; panel B, number of grids per side and panel C, 
the box side riox- In the last panel (D), we show the results 
of changing the guess value for the power spectrum Pw{k) 
in the weight function using equation (9). In panel A it can 
be appreciated the effects of three different weight assign- 
ment schemes: nearest grid point (NGP, solid line), cloud in 
cell (CIC, dotted line) and triangular shaped cloud (TSC, 
dashed line) . These weight assignments produce similar win- 
dow functions although the discreteness smearing effects in 
the high-order scheme are nontrivial (Jing & Valdarnini, 
1993). In panel B we show {\W{k)\/Nr)^ using different val- 
ues of A = 256 (solid Unc), 128 (dotted Une) and 64 (dashed 
line). As it can be seen, at large wave- numbers the increase 
of N produces a smoother window function. In panel C we 
display the results for different box sides rbox = 2, 3 or 4 
Teat where rcat is the minimum box side that contains the 
total cluster catalogue {rcat = 1037.1/i~^Mpc for XBACS 
and rcat = 1384.9/i~^Mpc for BCS). As it can be seen in 
this panel, they are indistinguishable. In panel D we show 
the dependence of the window function on the guess value, 
Pn,{k) =150000 (solid line), 500000 (dotted line) and 50000 
/i~^Mpc^ (dashed line). The window function amplitude is a 
monotonally decreasing function of Pw(k). This dependence 
is relevant for the power spectrum determination of BCS and 
almost negligible for XBACS (see panel D in Figure 5). We 
recall that the FKP method is valid for wave-numbers where 
the window function is a factor 2 lower than its value at the 
origin (Tegmark et al. 1998). In our case this corresponds 
approximately to A: > 0.01/iMpc~^. 

In Figure 5 we show the estimated redshift space power 
spectrum where the 4 panels correspond to those of Fig- 
ure 4. To make this plot clearer the logarithmic amplitude 
of BCS power spectrum was shifted -0.5. Panel A shows 
a strong suppression of power at high frequencies which is 
more significant than in the window function plotted in Fig- 
ure 4A. At low wave-number the agreement between the 
three schemes (NGP, CIC and TSC) is remarkably good. 
The increase of the spatial resolution using a larger number 
of cells per side produces a smoother power spectrum which 
can be appreciated in panel B. We have changed the guess 
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Figure 6. The upper panel shows the power spectrum for the 
complete XBACS (open circles) and the power spectrum for the 
XBACS clusters with z < 0.09 (filled circles). The lower panel 
shows the power spectrum for the complete BCS (open circles) 
and the power spectrum for the BCS cluster with z < 0.11 (filled 
circles) . 
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Figure 7. Comparison of the redshift space power spectrum of 
X-ray clusters with other samples obtained by different authors 
as indicated in the labels. 
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value P^{k) in the range 50000 /i^^Mpc^ - 
finding that the results are not strongly dependent on the 
particular choice of this value. We also find that the re- 
sults do not change substantially for box sides rhox='2, 3 
and 4 r^at- Accordingly, wc have used the NGP weighting 
assignment with parameters N = 256, rbox = 2rcat and 
Pw{k) = ISOOOO/i^^Mpc^. We have also tested the depen- 
dence of the window function and power spectrum results 
on the method used (equation (3) or equation (4)) to esti- 
mate the cluster rmmbcr density n(r) finding no significant 
differences in the results. In what follows we adopt the FKP 
fit (equation (3)). 

In order to test the presence of systematics effects that 
could artificially increase the amplitude of the power spec- 
trum at large scales {k < 0.03hMpc~^) (but see Miller & Ba- 
tuski 2000), we have restricted our samples to smaller sub- 
samples. Following Plionis & Kolokotronis (1998) we have 
generated two subsamples introducing a cut at z < 0.09 for 
XBACS and z < 0.11 ECS. In Figure 6 we show the corre- 
sponding power spectrum of each subsample compared with 
the power spectrum of the complete sample (Upper panel 
for XBACS and lower panel for BCS). From this compari- 
son wc observe a decreasing amplitude of the power spec- 
trum in both subsamples mainly in large scales. Since we 
know that in large scales the fluctuation signal have a very 
low level, even small systematic errors inherent to the survey 
construction can cause and overestimate of the power spec- 
trum in that scales. Consequently, wc adopt this subsamples 
like a more accurate estimator of the power spectrum for the 
XBACS and BCS. 

In Figure 7 we plot the redshift space power spectrum 
of XBACS (filled circles) and BCS (filled triangles). We have 
also computed the power spectrum of the recently released 
EBCS (filled squares) for z < 0.11. We compare our esti- 
mates with results obtained by other authors for different 
(optical and X-ray) cluster samples. The power spectrum 
obtained for XBACS has a higher amplitude than any other 
sample mainly in the range k > 0.05hMpc^^ . 

Since in Figure 7 we find different amplitude between 
XBACS and BCS we have tested if this discrepancy could 
be explained by the fact that both samples have different 
limiting fiux. Moscardini et al. (2000) show that for a given 
sample varying the limiting flux change the amplitude of the 
power spectrum but not the shape. In order to investigate 
this point, we have taken a subsample of the BCS cata- 
logue with the same redshift and fiux limit that XBACS 
and computed the power spectrum of this subsample. The 
obtained power spectrum for this subsample is very simi- 
lar to the power spectrum for the complete BCS. From the 
comparison wc conclude that the small variation in the lim- 
iting fiux do not explain the different amplitude observed 
in Figure 7 between XBACS and BCS. Moreover, EbeUng 
et al. (1998) suggest that different fiuxes could be assigned 
to a given cluster due to for example to different correction 
factors. This kind of differences can cause clusters to have 
the flux limit of 5.0xl0~^^erg cm^^/s when the BCS fluxes 
are used, while they remain just below the flux limit in the 
XBACS. Another possibility to explain the different am- 
plitudes is that XBACS could be richer and probably more 
massive than BCS that includes Abell, Zwicky and X-ray se- 
lected clusters with a poor optical counter-part. The power 
spectrum of XBACS is in good agreement with the power 
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Figure 8. Comparison of the redshift space power spectrum of 
XBACS (filled circles) and BCS (filled triangles) with theoretical 
models computed by Moscardini et al. (2000) as indicated in the 
labels. The error bars are computed using equation (13) for k < 
0.1/iMpc~^ and bootstrapping technique for k > 0.1/iMpc~^. 



spectrum of the REFLEX (Schuecker et al. 2001) cluster sur- 
vey for k ^ 0.03. We do not flnd strong evidence of a peak 
at wave-rmmbcrs ~ 0.05/iMpc^^, consistent with the results 
found by Miller & Batuski (2000). We find a fiattening of 
the power spectrum at these wave-numbers, smoother for 
XBACS than for BCS. The power spectrum of EBCS has a 
shape comparable to the original BCS except in large scales. 
In that scales, the BCS power spectrum have a higher am- 
plitude. 

In Figure 8 wc present the power spectrum error bars of 
XBACS and BCS computed using the N-body simulations 
(section 3.3). This errors have comparable size to errors ob- 
tained using equation 13 even they are internal errors that 
do not take into account any source of systematic errors or 
cosmic variance. We also plot the power spectrum predic- 
tions given by Moscardini et al. (2000) for the two sam- 
ples in two different cosmological models: Open Cold Dark 
Matter (OCDM) and Lambda Cold Dark Matter (LCDM). 
Both models have a dimcnsionlcss matter density param- 
eter flm = 0.3 but a dimensionless cosmological constant 

= for OCDM and Qa = 0.7 for LCDM. The normal- 
ization of the models is in agreement with observed cluster 
abundances (root mean square mass fluctuation on a sphere 
of radius 8/i"^Mpc, as = 0.87 for OCDM and as = 0.93 for 
LCDM). Moscardini et al. (2000) also present the SCDM, 
rCDM and TCDM models, but due to their signiflcantly 
lower amplitudes we do not included them in the compar- 
ison. From this plot it is possible to note that the obser- 
vational estimates for XBACS arc consistent in shape with 
the prediction of Moscardini et al. (2000) model for both 
OCDM and LCDM but with an amplitude ~ 1.4 higher 
that the models. That is probably an indication that these 
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Figure 9. In the left panel we plot the two point correlation 
function of XBACS (filled circles) applying equation (17) to the 
previously determined power spectrum. We also show the results 
of Abadi, Lambas & Muriel (1998) (open circles) and Borgiani, 
Plionis & Kolokotronis (1999) (open triangles). In the right panel 
we plot the two point correlation function of BCS using equation 
(16) (open circles) and (17) (filled circles). 

clusters that have a previous optical identification arc very 
massive. For BCS the agreement with theoretical predictions 
is better in the range k > 0.03hMpc~^ . As pointed out in 
Section 3 we can not take into account the ROSAT expo- 
sure time across the survey, i.e. we not use the sky coverage 
of the satellite because this is not available for any of this 
catalogues. We know that variation can affect the measured 
power spectrum especially on large scales. Consequently, this 
could bo the reason of a very high amplitude of BCS power 
spectrum at largo scales (k < O.OShMpc^^). The reality of 
the 'big bump' at A: = 0.015hMpc~^ could be analyzed in 
more detail when the survey selection functions of the sam- 
ples are taken into account properly if there were publically 
available. Moreover, the inhomogonous selection of clusters 
in the (E)BCS samples as well as the incompletenesses of 
XBACS and eBCS of 80 percent and below 75 percent, re- 
spectively, could introduce some artifical fluctuations. This 
might give also an explanation of the high fluctuation power 
detected on scales larger 400h~^ Mpc compared to LCDM or 
OCDM models. 

In order to test our results we derive an estimate of 
the two point correlation function from the power spectrum 
obtained previously. There are different ways to estimate the 
two point correlation function of a sample. One possibility 
is to derive it directly form the spatial distribution (Peacock 
1999): 

e(r) = {DD)/{DR) - 1 (20) 

where {DD) refers to the number of pair data-data sep- 
arated by distance r and {DR) refers to the same quan- 



tity for data-random pairs. This approach is the one applied 
by Abadi, Lambas & Muriel (1998) and Borgani, Plionis & 
Kolokotronis (1999) to XBACS. Other possibility is to es- 
timate 5 (r) using the relation between the power spectrum 
and the correlation function as Fourier transform pairs: 

^^r) = ^.jym^'-^dk. (21) 

In Figure 9 we have compared these two methods to 
estimate (,{r). For XBACS (left panel) we show the deter- 
mination of ^(r) using equation (16) given by Abadi, Lam- 
bas & Muriel (1998) (open circles) and Borgani, Plionis & 
Kolokotronis (1999) (open triangles). We also plot our do- 
termination of 5(r) using the estimation of the power spec- 
trum shown in Figure 6 though equation (21) (filled circles). 
This plot shows that there is a very good agreement between 
these three determinations. For BCS (right panel) we have 
computed ^(r) using both methods: equation (20) (open cir- 
cles) and equation (21) (filled circles). We find a reasonably 
agreement between the methods. Error bars of ^(r) are ob- 
tained from the power estimation uncertainties propagated 
through equation (21). 



5 CONCLUSIONS 

In this paper we apply the method proposed by FKP to 
compute the redshift space power spectrum of two flux lim- 
ited X-ray cluster samples: the XBACS and the BCS. At 
large wave- numbers 0.05/iMpc~'^ ^ 0.3/iMpc~^ we find 
a power law behavior P{k) oc /c" with n ~ —1.2 for XBACS 
and n ~ -1.0 for BCS. The power spectrum of XBACS 
is consistent with that of the REFLEX cluster sample. It 
should be remarked that the power spectrum amplitude of 
XBACS and REFLEX are significantly higher than the one 
derived for optical samples, which are quite consistent with 
our BCS determination. We do not detect the presence of 
any strong peak near 0.05/iMpc^^ supporting the idea that 
the presence of such a feature could be associated with bi- 
ases in the samples (MUler & Batuski 2000). The OCDM and 
LCDM models of Moscardini et al. (2000) produce a good 
fit for BCS over a long range of wave-numbers, but under- 
estimates tfie amplitude for XBACS. The shape of any of 
these models is also inconsistent with the shape obtained 
for XBACS. The correlation function of XBACS and BCS 
derived from the power spectrum agree with direct deter- 
minations (Abadi, Lambas & Muriel 1998 and Plionis & 
Kolokotronis 1998, this work). 
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